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Abstract 

In this paper, we will give a sufficient condition for a non-negative random variable X to 
Oh ' be heavy tailed by investigating the Laplace-Stieltjes transform of the probability distribution 

(~| . function. We focus on the relation between the singularity at the real point of the axis of 

convergence and the asymptotic decay of the tail probability. Our theorem is a kind of Tauberian 
theorems. 
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Q ■ 1 Introduction 
O 

ON 

Q ■ We consider the asymptotic decay of the tail probability P{X > x) of a heavy tailed random variable 
On : 

, X. A random variable X is said to be heavy tailed if 

V logP(X>x) ^ 

k>( , 7/ = hm / < 0. (1) 

x^QO log X 

7] is called the decay rate of P{X > x). 

Let F{x) be the probability distribution function (pdf) of X, i.e., F{x) = P{X < x). For 
example, if F{x) = 1 — x>l, then P{X > x) = hence X is heavy tailed with decay rate 
7] = —1. 

In this paper, we will give a sufficient condition for X to be heavy tailed by analytic properties 
of the Laplace-Stieltjes (LS) transform of F{x). The LS transform of F(x) is defined by 

POO 

^(s) = / e— (2) 
Jo 

In general, for a function R{x), which is of bounded variation in the interval < x < c for any 
positive c, the LS transform ^'(s) = e~^^dR{x), s = a + ir, is defined. If ^'(s) converges for 
fj > (To and diverges for a < uo, then do is said to be the abscissa of convergence of ^'(s). The line 
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5Rs = do is called the axis of convergence. In the case (Tq = 0, some local information of ^'(s) at 
s = provides the asymptotic behavior of R{x) as x ^ oo. Such a proposition is called a Tauberian 
theorem. The following is one of the Tauberian theorems. 

Theorem (Widder [14J, p. 192, Theorem 4.3) Let R{x) be a non- decreasing function and the abscissa 
of convergence of ^{s) be ao = 0. If for constants r < and A 

lim \'^{s)s-'' - A\ = 0, (3) 

then 

lim \R{x)r{r + l)x'' - A\ = 0, (4) 

where T denotes the gamma function. 

In [S] , [H] , [12] ! we studied the asymptotic decay of a light tailed random variable. A random 
variable X is said to be light tailed if the tail probability P{X > x) decays exponentially, i.e., 

lim -logP(X >x) <0. (5) 

We obtained in [12j the following theorem which gives a sufficient condition for a light tailed random 
variable. 

Theorem (Nakagawa [T2]) For a non-negative random variable X with probability distribution func- 
tion F{x), let ip{s) = Jq°° e~^^dF{x) be the Laplace- Stieltjes transform of F(x) and ctq be the abscissa 
of convergence of (p[s). We assume —oo < ctq < 0. If s = ao is a pole of if{s), then we have 

1 

lim - log P{X > x) = do. (6) 

X— >oo X 

It is known that for a monotonic R{x), s = ctq is a singularity of ^{s) (see Widder [T3], p. 58, 
Theorem 5b). Since pdf F{x) is monotonic increasing, ctq is a singularity of ip{s). If F{x) is the pdf 
of a heavy tailed random variable, then the abscissa of convergence of (p{s) is necessarily ctq = 
(see Widder [Uj, p. 40, Theorem 2.2b). Since 93(0) = dF{x) = 1, ctq is not a pole, but other 
type of singularity. 

In the research of the asymptotic decay of a tail probability, we would like to construct a general 
theory such that a local analytic information of ip{s) at s = fio tells the asymptotic evaluation of the 
tail probability. In a light tailed case [9].|llj.[T2]. we applied to this problem Ikehara's Tauberian 
theorem [5J,[7J and its extension Graham- Vaaler's Tauberian theorem [3], [7]. Ikehara's theorem 
assumes a global analytic property of 92(5), that is, s = ctq is a pole and there exist no other 
singularities on the axis of convergence = fio. While, Graham-Vaarler's theorem only assumes 
s = (To is a pole, which yields weaker assertion than Ikehara's theorem, however, it is enough for 
our purpose to investigate the asymptotic decay of a tail probability. In this paper, we will apply 
Graham- Vaaler's Tauberian theorem to the decay of heavy tailed random variable. 
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The purpose of our study is to apply our theory to the performance evaluation of the packet 
network. According to the research results for the Internet packet stream, it is reported that many 
characteristics are approximated by heavy tailed random variables [Tj. Contrary to the light tailed 
case, if, for example, the packet length is heavy tailed, the network performance, such as packet loss 
probability or end-to-end packet delay becomes worse. So, it is very important to investigate the 
tail probability of heavy tailed random variable from the view point of network engineering. When 
we apply the queueing theory to network engineering, even if we do not obtain the tail probability 
of packet length explicitly, we may obtain its Laplace transform by algebraic manipulation as 
Pollaczek-Khinchin formula [6]. Then, we know the singularity of the Laplace transform and the 
asymptotic decay can be investigated by our theory. 

Throughout this paper, we will use the following symbols. N, N^, M, C, denote the set of 
natural numbers, positive natural numbers, real numbers, complex numbers, and further, SR, £, T 
denote the real part of a complex number, Laplace transform and Fourier transform, respectively. 

2 Examples of Heavy Tailed Random Variable 

Now, we look at some examples of heavy tailed random variables and their LS transforms. 

2.1 Continuous Random Variable I 

Let X be a random variable with pdf F{x) = 1 — 1/x, x>l. The decay rate is r/ = —1. The LS 
transform 93(5) of F{x) is represented in a neighborhood of s = as 

/oo 
e— (7) 

= slogs + /3(s), (8) 
where /?(s) is analytic in a neighborhood of s = (see Lemma [T]). 

2.2 Continuous Random Variable II 

For a pdf F{x) = 1 — 1/ ^/x^ x > 1, we have the decay rate r] = —1/2 and 

<^(^) = + (9) 

where /3(s) is analytic in a neighborhood of s = (see Lemma [2]). 

2.3 Discrete Random Variable 

Let X be a discrete random variable with probability distribution p = {pn)n€N] 

= - — L_ - - — L_, n G N, r e N+. (10) 
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The tail probability of X is P{X > x) 
generating function (pgf) p{z) of p is 



(n + 2y 



and the decay rate is ry = — r. The probability 



p{z) 



z - 1 -e- 1 

n=l 



Substituting z = e we have 



^is)=pie-') (12) 
= a(s)sMogs + /3(s), (13) 

where a(s) and /3(s) are analytic in a neighborhood of s = with a(0) / (see Lemma [3]). 

2.4 Stationary Distribution of M/G/1 Type Markov Chain 

We consider the tail probability of the stationary distribution of an M/G/1 type Markov chain 
[2j. Let X be a heavy tailed discrete random variable with probability distribution b = (6n)neN- 
Denote by B{z) the pgf of b. Consider an M/G/1 type Markov chain with probability transition 
matrix 



P 



( bo bi 62 ^3 

ao ai 02 as 

ao ai a2 

ao ai 

V : : ; ; 



\ 



(14) 



where a = (a„)„gN is a probability distribution with pgf A{z). Suppose there exists the stationary 
distribution vr = {■7Tn)neN of P. Let 7r(z) = X^^o '^nz'"' be the pgf of vr, then by the Pollaczek- 
Khinchin formula p].[3].[T0]. 

^0 {zB{z) - A{z)) 



tt{z) 



(15) 



z-A{z) 

Substituting z = e^^ into (jlSp . we have 99(5) = 7r(e^*). If the singularity of tt{z) at z = 1 comes 
from the singularity of B{z), we know the singularity of f{s) at s = comes from B{e~^). So, it 
is expected that vr should be heavy tailed. 

3 Main Theorem 

From above examples, we expect that the following theorems hold. These are main theorems in 
this paper. 

Theorem 1 Let X be a non-negative random variable with probability distribution function F{x), 
and 99(5) be the Laplace- Stieltjes transform of F{x). Assume the abscissa of convergence of (p{s) is 
(Tq = 0, and ip{s) is represented in a neighborhood of s = as 



(p{s) = a{s)s''logs P{s), r e N+, 



(16) 
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where a{s) and (3[s) are analytic with a(0) 7^ 0. Then, we have 

lim = (17, 

x^co log X 

Theorem 2 Under the same notation as in TheoremUl if ctq = ^ o-i^d f{s) is represented in a 
neighborhood of s = as 

ip{s) = a(s)s^ + /?(s), r > 0, r ^ N, (18) 
where a{s) and /3(s) are analytic with a{0) / 0. Then, we have 

x^oo log X 

The following is a related work to Theorem 2. 
Theorem (Korevaar [7\, p. 194, Theorem 8.2) Let S{x) vanish for x < and be locally integrable, 
positive and non-increasing for x > 0. Let 99(5) be the Laplace- Stieltjes transform of S{x). Further, 
let p{x) = x~'''l{x) with < r < 1 and slowly varying l{x). Then, for a constant A, 

Six) — S(oo) , , , 

lim , . = ^ (20) 

i'-»oo pi^x) 



if and only if 



^,^m^h:IM = ATil-r). (21) 

p[x) 



x^oo 



4 Preliminary Lemmas 

We will prepare some lemmas for the proof of our main theorems. First, let us define the step 
function Ai(t) as 

Lemma 1 For r £ N^, we have 

^{s)^c(^j^A^{t)^ (23) 
-e-'^dt (24) 



00 1 

St, 



where j3{s) is analytic in a neighborhood 0/5 = 0. 



^ ^^\^logs + (3{s), (25) 
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Proof The (r + l)th derivative of (p{s) is 



POO 

= (-If+W e-''dt (26) 



1 

-irM'^+^r (27) 



Since (e — l)/s is analytic, we have, by successive integrations, the desired result. □ 

Lemma 2 For r > 0, r ^'H, let = \r\ he the maximum integer among integers smaller than r, 
and let f = r — rQ. Then we have 



(^(s)^£(^Ai(t)) 
(_l)ro+V 



r(r + 1) sin vrr 
where P{s) is analytic in a neighborhood of s = 0. 



(28) 

s' + Pis), (29) 



Proof By a formula of Laplace transform (see [8], p. 287), 

/>oo 

<^('-o+i)(s) = (_i)-o+i / t-^e-''dt (30) 



1 

/•OD 

= (-1)'°+^ \ - I \ t-'e-'*dt (31) 
= (_l)-o+ir(l - f)/-i + /3(s), (32) 

where /3(s) is analytic in a neighborhood of s = 0. By successive integrations, 

^{s) = {-ly^+^ni -r)-_.^-...- 3^/+^o + ^(,) (33) 

r r + 1 r + ro 

(_l)ro+l7r 

:s" + /3(s). (34) 



r(r + 1) sin vrr 

In (fM]l . we applied the formulas; r(2; + 1) = 2;r(z) and r(z)r(l — z) = vr/ sinvrz. □ 

Lemma 3 Let X he a discrete random variable with probability distribution p = (pn)neN; 

Pn = 7 — - 7 — n G N, r e N+, (35) 
(n + 1)'- (n + 2)^' ' ^ ' 

and let p{z) = Y2'n=oPn^^ Pdf (^f P- Then, f{s) = p{e~^) is represented in a neighborhood of 

s = as 

V9(s) = a(s)sMogs + /3(s), (36) 
where a{s) and /?(s) are analytic with a(0) / 0. 
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Proof For r > 2, by calculation and Riemann's formula (see Widder |14j . p. 232), 

n=l 

1 z-l f°° zf-^ ,1 

Then, 

By the change of variable u = e*"*"* — 1, we have for sufficiently small s > 0, 

e^(l-e^) {log(l + -u) - s}^-i ^ ^ 

cp(s) = / ^ — ^— ! dtt + regular term (40) 

Je^-l U{U + 1) 

[ (-- {{-sy-' + (r - l){-sy-\ +...}du + r.t. (41) 

\r-l ' 1 



e"(l - 




(r - 


-1)! 


e"(l - 




(r - 


-1)! 


e"(l - 




(r - 


-1)! 



where 



i-sy-' - + r.t. (42) 
a(s)snogs + /3(s), (43) 



«(-) = V^^'^' = ^ ' ^''^ 

and /?(s) is some analytic function. 
For r = 1, we have 

Piz) = ^^iog{l-z) + ^. (45) 
Similar argument leads to the desired result. □ 
Lemma 4 For a pdf F{x), let 

POO 

ip(s) = / e-'^'dF^x) (46) 
have the abscissa of convergence ctq = 0. If, in a neighborhood of s = 0, 

ip{s) = a{s)s'' log s + /3(s), r G N+, (47) 

where a{s), /?(s) are analytic with a(0) ^ 0, then, 

r > 0, if r is odd, 
< U, II r IS even. 



Proof From (IM 



> 0, II r IS even. 
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while from (j47p by calculation 

(f'^'^'' (s) = r\a{s) log s + 0{s), 
where 9{s) is a function of s with |0(O+)| < oo. Thus, by ([SD]) 

(^M(0+) = r!a(0) x (-00) + 9{0+). 
Comparing (09]), ([5T]) . we have the result. 



(50) 

(51) 
□ 



Lemma 5 Under the same notation as in Lemmas if fi^) has the abscissa of convergence 

(To = and 



ip{s) = a{s)s'' + P{s), r > 0, r ^ N, 



with a(0) 7^ 0, t/ien, 



Proof From (IM 



while from (|52 



a(0) 



^(ro+i)(o+) 



> 0, if ro = [rj is odd, 
< 0, if ro is even. 



> 0, if ro is odd, 
< 0, if ro is even. 



ro+l 



(r-o+1) 



k=0 



+ ^(r-0+l)(, 



Since 



0, for /c = 0, 1, . . . , ro 

s=o+ 1 +00, for A; = ro + 1, 



(52) 
(53) 
(54) 

(55) 

(56) 

(57) 
□ 



we have from 

^(r-o+i)(o+) = a(0) X (+00) +/?(^o+i)(0)- 
Comparing (j^Tj) . we have the result. 
4.1 First Several Terms of a{s) and P(s) 

We will need later, in the proof of main theorems, to make Laplace transforms which have the same 
first several terms as those in the Taylor expansion of a{s) and /?(s), respectively. 
First, consider the following case 



ip{s) = a{s)s'' log s + /3(s), r G N"' 



(58) 



with expansions a{s) = Xl^o'^"-^"' '^(0) 0' ^^'^ (^(^) ~ YlnLo f^nz"' ■ 

We will make functions g*{t) and h*{t) such that their Laplace transforms G*{s) = C{g*{t)) 
and -&*(«) = C{h*{t)) satisfy the following (i) and (ii) for L G N+. 
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L-1 



(i) G*{s) = a*{s)s'^logs + (3{s), where a*{s) = an^" and (3{s) is analytic in a neighborhood 

n=0 

oo 

of s = 0. Let (3{s) = be the expansion at s = 0. 



n=0 



(ii) ii"*(s) = ^(/?n - /3n)s"+ higher order terms. 



n=0 



The above (i) and (ii) mean that a*(s) is equal to the sum of the first L terms of a(s), and the 
sum of the first L terms of H*{s) + (3{s) is equal to that of (3{s). 
Function g*{t) 



Define 



L-l 



k=0 
\r+k+l 



9k = (-ir+'=+^(r + ky.ak, k = 0,1,..., L-l. 



(59) 
(60) 



where Ai(t) was defined in ()22p . By Lemma [H we see that g*{t) satisfies (i). The first coefficient 
go is positive by Lemma HI i.e., 



go = (-l)^+V!ao > 0. 



(61) 



Function h*{t) 

Let hk{t) = ke~''\ t>0, = 1, 2, . . . , and Hk{s) = C{hk{t)). We have Hk{s) = k/{s+k), > 
—k, and the expansion 

oo 

/ fi \ TJ. 

(62) 



n=Q 



The sum of the first L terms of (|62p is represented as 

L-l\ 



1 



1 



L-1\T 



[l,s,...,s^-') 



(63) 



where denotes the transposition of vector. Let V be the L x L matrix; 

/ 1 -1 ... (-1)^-1 \ 

L-l 



V 



1 



1 

2 

1 
L 



L-l 



(64) 



We will make a desired function by a linear combination of hk{t), k = 1,2, ... ,L. Let h{t) = 
'nik=idkhk{t), t > 0, and write d = {di,d2, . . . ,dL), s = {1, s, . . . , s^"^). Further, write (3 = 
(/3o, . . . , /3l-i), /3 = (/3o) • • • ) /3l-i)- The sum of the first L terms of H{s) = C{h{t)) is 
cLYs^ , then we must solve the equation 



dys^= (/3-/3)s^ 



(65) 
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Since det F 7^ (Vandermonde matrix), we have d = (P — (3)V ^. We write this solution as 
d= {di,d2, ■ ■ ■ , di), then 

L 

h*{t) = ^dkhk{t) (66) 

k=l 

is a desired function, i.e., H*{s) = C{h*{t)) satisfies (ii). 
Summarizing above. 

Lemma 6 Let (/?(s) be the LS transform of a pdf and the abscissa of convergence be Uq = 0. If 

(^(s) = a(s)sMogs + /3(s), r G N+, (67) 

where a{s), /5(s) are analytic in a neighborhood of s = with a(0) 7^ 0, then g*{t) in l\59t\ and 
h*{t) in satisfy (i) and (ii). 

Similarly, in the case 

(^(s) = a(s)s'" + /3(s), r > 0, r ^ N, (68) 

with a{s) = XlJ^Lo '^n'^"' "^0 7^ 0) /5(s) = Z]^o'^"'^"' ^^^^ make functions (7*(t), h*{t) such that 
their Laplace transforms G*{s), H*{s) satisfy the following (iii) and (iv) for L G N^. 

L-l 

(iii) G*{s) = a*{s)s'^ + (5{s), where a*{s) = a„s" and f3{s) is analytic in a neighborhood of 

n=0 

00 

s = 0. Let (3{s) = /3ns" be the expansion at s = 0. 

n=0 

L-l 

(iv) H*{s) = ^(/3n - Pn)s"'+ higher order terms. 

n=0 

In this case, 

L-l 

5*W = E^^^iW' (69) 

fc=0 

5fc = {-lyo+k+i^J^^^ + l)ak, k = 0,1,..., L-l. (70) 
vr 

satisfies (iii). The first coefficient go is positive by Lemma [5l i.e., 

50 = (-l)^o+i^^r(r + l)ao > 0. (71) 

TT 

The same /i*(t) as ()66p satisfies (iv). Thus, we have 

Lemma 7 Lei (p{s) be the LS transform of a pdf and the abscissa of convergence be o"o = 0. If 

ip{s) = a(s)s^ + P{s), r > 0, r ^ N, (72) 
then g*{t) in <\69^ and h*(t) in l\66\i satisfy (iii) and (iv). 
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4.2 Majorant and Minorant Functions 

For the evaluation of the tail probability P{X > x) from above and below, we need to use majorant 
and minorant functions for an exponential function (see Korevaar [7J, p. 132, Graham- Vaaler [3]). If 
two functions /i, /2 satisfy fi{t) > f2{t), t £ M, then /i is said to be a majorant for /2, and /2 is 
a minorant for fi. 

For LO > 0, we will define a majorant M^(t) and a minorant ml^^t) for 

f e~'^* t > 

^"<" - { 0, ' * ; 0. 

Define (see Korevaar [7J, p. 132) 

M^{t)=(^-^YQUt), teR, (74) 



vr 



n=0 n=l ^ ' 

and 

mi(t) = Mi(t)- (^y, tGM. (76) 

Moreover, for L G define 

M^{t) = {Ml{t)f and m^(t) = {rri{t)f . (77) 
For fj > 0, 5 > 0, write uj = 2T:a/5, then define 

M.Vt).M.^(|)^M4./.(|), (78) 



= "if; (^) = [^]- (79) 



(2vr) ■ 



Lemma 8 (Korevaar [7]) For any L G N"*", o" > 0, (5 > 0, 

Mi;^, m^_5GLi(M)nL2(M). (80) 

For A > 0, an entire function f{z) of a complex variable 2; = x + is of exponential type A if 

\f{z)\ < Cexp(A|z|), zeC, C>0. (81) 

A real function f{x) is of type A if f{x) is the restriction to M of an entire function of exponential 
type A. 

Lemma 9 (Korevaar [7], Nakagawa [T2]) and are of type L6. 
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Lemma 10 (Korevaar Graham- Vaaler [3]) For L G N^, 



and for odd L £ 



(82) 



(83) 



Proof If L = 1, the result fohows Korevaar |7], p. 129, Proposition 5.2. The odd power preserves 
the order of real numbers. □ 
From Lemma [51 we can define the Fourier transforms M^^ = !F{M^^) and = J-{m^g), 



where the Fourier transform is defined as 



(84) 



Then, from Lemma E] and the Paley- Wiener theorem, we have 



Lemma 11 (Rudin [13], Korevaar [7J) For any L £ N^, 



supp{M^g) C [-LS,L6] and supp(m^5) C [-L(5, L(5], 



(85) 



where supp denotes the support of a function. 



4.3 Calculation of M^^ and m^^ 



It is not difficult to calculate the Fourier transforms = J^(M^) and m^J^ = J^{mlj). 
Define 



sinvrt 



92 (t) 



sin'^ vrt 



irt J Tit 
and write qi = J^{qi), q2 = •^('72) • By calculation, we have 

T 



t G 



(86) 



1 + — , -27r<T<0, 

1-7^, 0<r<2^, 
zvr 



(87) 



I 0, 



and 



(l2{r) 



2' 



otherwise. 



-27r < r < 0, 



< r < 27r, 



0, otherwise. 
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Lemma 12 We have 

1 _ e-('^+«T) ' ~ n I 1 - e-('^+«^) ^ 1 



and 



Proof See Appendix 1X1 

Next, we will calculate = J^{M^), rh^ = J^{m^) and then calculate Unit, 
lim^^o+ rh^{T), for r / 0. 



Let us define 



(\ 2 9 
sm vrt \ io sm 7rt 
vrt y IT 

U) 

— ( 

vr 



- -92(t), tern, 



oj sin^ vrt 
vr vrt 

-q2{t), t £ 

TT 



and Uuj = T{u^), Vuj = 
Lemma 13 We have 



where * denotes the convolution operation and *l denotes the l-fold convolution. 
Proof See Appendix [Bl 
Lemma 14 For r 7^ 0, 



L -I / -r \ / — iV \ ^ 



Proof (ETI) and (EHl) follow 



= E ^ (^) (^) «V) . a-V). 



u 

' it 



v:'~\r) = (^) " #-'(r). 
By the change of variables, we have 
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Lemma 15 For r 7^ 0, 

0-^0+ 



27r 1 1 /L\ f e-^2W5 V /27rr\ , /27rr\ 



5 Proof of Theorem [T] 
5.1 Upper Bound for P(X > x) 

First, we will evaluate P{X > x) from above by using the majorant function M^^. 
Let L G N+ with L>r. For arbitrary cri > 0, (72 > 0, (5 > 0, 



/•oo 

= / EL,,(t-x)e-'^i*(iF(t) (104) 
Jo 

POO 

< / Mi; sit - x)e-"i*o!F(i), x > 0, (105) 

JO 

where the last inequality holds by (182p in Lemma [10] (see also Korevaar [7] , Nakagawa [12j ) . By 
Lemma [TT\ M^^ ^(t — x) is represented by the inverse Fourier transform of M^^ 5 — ("^<t2 s) ^ 

KAi - ^) = ^ y^^^Mi;_,(-T)e*(--*)-dr. (106) 



(107) 



Substituting (1106P into (jlOSp . we have by Fubini's theorem 

/"CXD /"OO -I /"Z/^ 

/ M^.it-x)e~''''dF{t)= — .(-r)e^("-*)^e"'^i*(irdF(t) 

= — / M^, 5(-r)e"^dr / e-('"i+^^)*dF(t) (108) 
2tt J_ls Jo 

= - y Mi;,5(-r)e-Xc7i + ir)dT. (109) 
Now, for the representation of v'(-5), i-e., 

V9(s) = a(s)s'" logs + /3(s), r G N+, (110) 

we define 

rit)=g*{t)+h*it), t>0, (111) 

where g*{t) and /i*(t) were defined in Lemma[6l Let f*{s) = £{f*{t)), then f*{t) and V'*(s) have 
the following properties (a), (b), (c) and (c'). 

(a) f*{t) > for all sufficiently large t. 
This is because g*{0) = an > bv (|UT]) . 
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Define ^(s) = ip{s) — ^p*{s), then by Lemma [H 

^(s) = ipis)-iG*is) + H*{s)) (112) 
( / oo \ oo ^ 

aL+nS"" s' log s + ^ (pL+n - Pl+n) > , s = a + ir. (113) 



I, \ n=0 / n=0 J 

We see from ()113p . 

(b) i^(s) is continuous in the closed region {0 < o" < e, —Ld < r < L5} C C for sufficiently 
small e,6 > 0. 

In fact, because s = is an isolated singularity, we can take €,5 so small that the closed region 
{0 < o" < e, —L5 < T < L6} does not include any singularities of (p{s) and ^*{s) other than s = 0. 

Define M^^si^) = lim^^o+ ^^^(r), rh^^sir) = lim<,^o+ -L6 < t < L6, r / 0, for 6 
sufficiently small as in (b). Then, we have 

(c) supp(M(^5) C [—L6, L6] and Mj^^(r)^(zT) is r times piecewise differentiable with 
[M,^^simir)f^ eL'{[-L6, L6]) . 

(c') supp(mQ^) C [—L6, L6] and rfiQ g{T)(,{iT) is r times piecewise differentiable with 
{m(^^,{r)C{iT)f^ eL'([-L6, L6]). 
These (c) and (c') hold because of Lemma [T5] and (jllSp . 

Now, defining F*{t) by dF*{t) = f*{t)dt, we have in a similar way as (jlOOp 

oo 1 /'L<5 

M4,5(t - x)e-'^-'dF*{t) = — / Mi;,5(-r)e->*(ai + ir)dr. (114) 

-^^ J-LS 

Subtracting pi4p from pU9p . we have 

poo poo 

/ M^^4t-x)e-'^^'dF{t)= M^^4t-x)e-''^'dF*{t) (115) 

JO ^0 

+ - / Mi;,5(-r)e--e(^i + ^r)dT. (116) 

For sufficiently small 5 > 0, 



^(fr) = lim ^(cJi +ir), -L(5 < r < (117) 

(71^0+ 



is uniform convergence due to (b), hence 

0-1^0+ 2tt J_j 
From (fTn5]l . (fTTHIl . (fUH]) . for cji ^ 0+, we have 



1 r^'' 1 f'^^ 

hrn - / M^^4-T)aai + iT)e^^dT = — / Mi;,5(-r)e(ir)e--dr. (118) 
^1^0+ 27r 27r J.^^ 



/"OO /"OO 1 P-L6 

e^^^^^l e-'^'^-'dFit) < M^^,,(t-x)dF*(t) + -y ^^Mi;,,(-r)C(ir)e--dr. (119) 
By the estimation for M^{t) (see Korevaar [7J, p. 132), i.e., 



1 , , / sinvrA 
< M^{t) < ( , * < 0, 



(120) 
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we have 



sin 



5t/2 \ 



2L 



6t/2 J 



sin vrt 



t < 0, 



t > 0, 



(121) 



where u; = 2'iTa/5. Thus, it is easy to see that there exists a constant Ci > such that 

2L 

( ' 

MUt-x)<{ 



0<t <x-l, 
t> x-1. 



^5{x-t)/2 

Therefore, the first term of the right hand side of (I119p is evaluated as 

/•oo 

,(t - x)dF*{t) = / M^Jt - x)nt)dt 



L-l 



(122) 



M!;^^,{t-x) {g*{t) + h*{t))dt 



(123) 
(124) 



1 



6{x-t)/2j 



2L 



px—1 

fc=0 



1 



fc=0 
L 



-dt + C^Y.\gk\ j^dt (125) 



+ (^(^3^) e-^'dt + C,^m\l_^e-^^dt (126) 



< ©(x-^'^+i)) + + 0(x-2^) + 0(e-"), C2 > 



< C3 > 0, 



(127) 
(128) 



for all sufficiently large x, by virtue of Lemmas I16t 1171 in Appendix. Notice 50 > and L >r. 

Next, the second term of the right hand side of ()119p will be estimated. We have by (c) and 
integration by parts. 



L5 1 rL5 



hin - / M^^4-T)air)e^''^dr = - / Mo^,(-r)e(ir)e--dr (129) 
CT2^o+ 27r J_Ls 27r J_ls 



2'KX- J_LS 

o{x~^), X — > 00 



(r) 



(130) 
(131) 



due to Riemann-Lebesgue theorem. Then in f)119p for a2 0+, we have by (1128P and ()13ip . 

poo 

P{X >x)= dF{t) < -, C> 0, (132) 

J X -X 

for all sufficiently large x. 

5.2 Lower Bound for P{X > x) 



We will evaluate P{X > x) from below by using the minorant function m^^. 
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Let L G be an odd number with L > r. For arbitrary fii > 0, (72 > 0, (5 > 0, 



La2X 



Jx Jo 



> 



m^^^g{t - x)e-'''^dF{t), x>0 



In a similar way as from (jl03p to (jllOp . we have 
Due to the estimation for mlj{t) (see Korevaar [7j, p. 132), i.e., 



sin nt 



TTt 



< mi{t) < 0, 



t < 0, 



-Ult 



sin vrt 
nt 



< ml{t) < e-"^* t > 0, 



we have 



sin5t/2\ 
6t/2 ) 



2L 



< m^Jt) < 0, 



-LUt 



sm 



5t/2 \' 



6t/2 J 



-Lujt 



t < 0, 



t > 0, 



where uj = 2iTa/5. Thus, there exit constants C4, C5, Cg > such that 

2L 



1 



-i)/2 
C4, 



5{x-t)/2 



< t < X - 1, 
rc-l<t<a; + l, 

, t>x + l. 
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Therefore, the first term of the right hand side of (jl35p is evaluated as 



-C,J2\9>^\ r^T^idt (140) 
fc=o -^^-1 

/'OO -1 

+c,y: bfci / ^-^'^'^'^^/'^^*-^^7^^* (141) 

^-1 /•°° / 1 \^ 1 



A:=0 



fc=l 

^ rx+1 



/o V^(^-i)/2 

-C4^A:|4| / ^ e~''*dt (144) 

^ poo 

+ C5^k\dk\ / e-(2-W'5)(t-^)e-fctrft (145) 



fc=l 



fc=i 



/>oo 1 

> ©(x-^'+i)) + ©(x-^'+i)) + C7^5o / e-(2-i-2/5)(*-^)__dt (147) 
+ 0(x-("+i)) + 0(x-2i) + 0(e-^) + 0(e-^) + ©(e"^), (148) 



where C5 is a positive constant. Thus, 



00 



hin / m^^^sit - x)dF*{t) >^,Cr> 0, (149) 

for all sufficiently large x. 

Next, we will evaluate the second term of the right hand side of (jl35p . We have, by (c') and 
the integration by parts, 

hm — / m^^,5(-r)e(ir)e"-dr = — / mo^,5(-r)e(ir)e--dr (150) 

= o(x~'^), X 00, (152) 
by Riemann-Lebesgue theorem. Then, in (jl35p for a2 0+, we have from p49p and ()152p . 

P(X >x)= dF(t) > — , C > 0, (153) 

for all sufficiently large x. 

From (|132p and (|153p . the proof of Theorem [T] is completed. □ 
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6 Proof of Theorem [2] 

The same proof as that of Theorem [1] is appUcable by replacing g*{t), h*{t) in Lemma [5] with those 
in Lemma [6j □ 

7 Conclusion 

In this paper, we investigated the asymptotic decay of the tail probability of a heavy tailed random 
variable. We proved two theorems which give sufficient conditions for a random variable to be 
heavy tailed. Our theorems are based on the Tauberian theorems due to Graham and Vaarler. The 
central idea is the approximation of the exponential function by a major ant and minorant functions 
whose Fourier transforms have a finite support. 

Through the proof of the theorems, I think that some more general representation of the sin- 
gularity guarantees the random variable to be heavy tailed. 
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Appendix 



Lemma 16 For ni,n2 € N"*", ni > 2, 712 > 2, let n = min(ni,n2). Then, 

< Oix-""), X ^ 00. 



dt 



{x - t)"it"2 
Proof By the change of variable t = xu, 

'■^-1 dt 1 /■^-^/^ du 



1 (x-t)"ir2 a-m+na+i y^^^ 



1 f /-l/^^ /-l"!/-] rf^ 



3,ni+n2+l 1 J^^^ J (1 - n)"iu"2 

< < 2"! / h 2"2 ' 



1 (' „?12 — 1 _ nn2~l ^ni — 1 

2"! h 2"2 - 



x"i+"2+i I n2 - 1 ni 

0(x~"), X ^ oo. 



Lemma 17 For /c > and n G N, we /laue 

r-a;— 1 ^—kt 



/ 7 — dt = 0(x~"), a; ^ oo. 



Proof By the change of variable u = x — t, 



(x-t)^ A 

= — / —du — 



1 

• , X ^ oo, 

/cx" 



from L'Hopital's rule. 

A Proof of Lemma [12 



■nui 



, \ vrft — n) J I vr I V vrft — n) 

n=o \ ^ ^ ' ' / n=o \ V / 

oo / oo oo \ 

e-"('^+^^) - ^92 (r) ^ e-'^('^+^^) - ^ e" 

n=0 \n=0 n=0 / 



1 _ e-('^+«'^) vr — e~('^+*'^) 1 — e 

The result for is proved in a similar way. 
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B Proof of Lemma [13 



/ sin vrA 



2L 



ni=0 nL=0 



+ 



{t - ulY t-riL t 
y. ^_„^^ f / sin7r(^-ni) \ ^ _ u; sin^ 7r(t - m) _^ a; sin^vrt 

I V vr(t — ni) / IT vrfi — ni) vr vrt 

ni=0 rti=0 I ^ \ ' ^ ^ ' 



oo oo 



sin7r(t — n/,)^ a; siii'^ 7r(t — n^) ^ cj sin^ vrt 



7r(t — ul) J IT ii{t — riL) TT Tit 



711=0 ni=0 



Therefore, 



oo oo ^ 

ni=0 ni=0 ^ 



oo oo 



(2vr)i-i E ••• E 

^ ' ni=0 n£=0 



X ^ I (e-^"^i^n.(T)) * ... * (e-*"'=<^n^(r)) ^ * l):^-'(r) (154) 

'=0 l,fci,...,fc; : distinct 

^ oo oo 

(27r)^-i E ••• E 

^ ^ ni=0 n£=0 

X E ( E ^"'"'^^^ • • • ^"''^''^ I (^) * *-^"'(^) (155) 

'=0 \A:i,...,A;; : distinct 



i=0 

oo oo 



ni=o ni=o \fci,...,fc, : distinct 



L\ / 1 V f 1 ^ 



\ , '■*L—l I 



C(r)*C-'(r). 



In (ji54p . p55p . p56p . for / = 0, the (empty) sum on ki, ki is considered to be 1. 
Similarly, we have the result for m^{T). 
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